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Abstract 

For a particular choice of the smoothing kernel, it is shown that the system of 
partial differential equations governing the vortex-blob method corresponds to 
the averaged Euler equations. These latter equations have recently been de- 
rived by averaging the Euler equations over Lagrangian fluctuations of length 
scale a, and the same system is also encountered in the description of inviscid 
and incompressible flow of second-grade polymeric (non-Newtonian) fluids. 
While previous studies of this system have noted the suppression of nonlinear 
interaction between modes smaller than a, we show that the modification 
of the nonlinear advection term also acts to enhance the inverse-cascade of 
energy in two-dimensional turbulence and thereby affects scales of motion 
larger than a as well. This latter effect is reminiscent of the drag-reduction 
that occurs in a turbulent flow when a dilute polymer is added. 



The two-dimensional incompressible, Euler equations are 



d t uj + V • (uw) = 0, V-u = 0, w(* = 0) = w , (1) 

where u> = V x u is the vorticity, it is the spatial velocity vector field, t denotes time, and 
all the dependent variables depend on t and x = (xi,X2), the Cartesian coordinates in the 
plane. An inversion of the vorticity- velocity relation yields u = J K(x, y)uj(y)dy, where 
K = 'V ± G, G is the solution of —AG = 5, and V = (— d Xl , d X2 ).For fluid motion over 
the entire plane, K(x, y) = (2n)^ 1 V -1 log \x— y\. Let r) t denote the flow of uk = u(t, •), so 
that drj t /dt = Ut(r](t)). Because ttj is divergence-free, the flow map r\ t is an area-preserving 
transformation for each t. It follows that 

^ = JK(r 1t (x),r ]t (y)Mr 1t (y))dy 

= J K(rj t (x), rj t (y))u (y)dy, (2) 

where the last equality is a consequence of the pointwise conservation of vorticity along 
Lagrangian trajectories, u(r) t (x)) = ujq{x). Thus, the initial vorticity field completely de- 
termines the fluid motion. Choosing the initial vorticity to be a sum of N point vortices Si 
positioned at the points x^ in the plane with circulations Tj, co> = Y^=\ r^i, equation @ 
produces the classical point- vortex approximation to (fl|). This approximation is known to 
be highly unstable, as finite-time collapse of vortex centers may occur.@ 

Chorin's vortex blob methodi alleviates the instability of the point-vortex scheme by 
smoothing each delta function Si with a vortex blob \i a function that decays at infinity, 
and whose mass is mostly supported in a disc of diameter a. Thus, instead of using the 
integral kernel K(x, y), one uses the smoother kernel K a = 'V ± G a where G a is the solution 
of —AG a = x- The vortex-blob method then evolves the point-vortex initial data, which we 
shall now call qo, by the ordinary differential equation 



dt 



J IT{^{x),^{y))q G {y)dy. (3) 



Henceforth, to keep the notation concise, we will drop the superscript a when there is no 
ambiguity. 

When the vortex-blob G a is the modified Bessel function of the second kind K , it is the 
fundamental solution of the operator (1 — a 2 A) in the plane, and the vorticity q is related to 
the smoothed velocity vector field u by q = (1 — a 2 A) V x u. Thus, Chorin's vortex method 
for this choice of smoothing is given by the partial differential equation 

d t q + V ■(uq) = 0, V-tt=0, q{t = 0) = q , (4) 

The system of equations (f|) is also known as the two-dimensional isotropic averaged 
Euler equations, and are derived by averaging over Lagrangian fluctuations of order a about 
the macroscopic flow field.B When the constant a > is interpreted as a material 
parameter which measures the elastic response of the fluid due to polymerization instead of 
as a spatial length scale, then (f|) are also exactly the equations that govern the inviscid flow 
of a second-grade non-Newtonian fluid.0 According to Noll's theory of simple materials, 
(H) are obtained from the unique constitutive law that satisfies material frame-indifference 



2 



and observer objectivity. Consequently, the vortex method with the Bessel function Kq 
smoothing naturally inherits these characteristics.!!! 

We find the connections between averaging Euler equations over Lagrangian fluctuations, 
a constitutive theory for polymeric fluids, and a classical numerical algorithm to be quite 
intriguing and suspect that these equations will be important from a modeling standpoint. 
However, most previous studies of the averaged Euler equations have been of a mathematical 
nature, and we are aware of only a few cases where this system has been used as a (dynamic) 
modeling tool: Chen et al.i used a viscous version of the three-dimensional averaged Euler 
equations to simulate isotropic turbulence and found that they could reproduce large-scale 
features without fully resolving the flow. Nadigai considered the inviscid two-dimensional 
form of the averaged Euler equations and demonstrated that for suitably chosen values of a, 
the large-scale spectral- scalings of the Euler equations could be preserved while achieving a 
faster spectral decay at the smaller scales. Finally, Nadiga & Marg olinS used an extension of 
the two-dimensional averaged Euler equations in a geophysical context to model the effects 
of mesoscale eddies on mean flow. 

Before we go on to consider numerical simulation of the averaged Euler system, we 
wish to point out that there is also a beautiful geometric structure to (Q) which follows 
the framework developed by ArnoldB and Ebin and Marsden.0 While the details of this 
particular issue are far outside the scope of this article, it is, nevertheless, worthwhile to 
state the result. Arnold showed that the appropriate configuration space for a perfect 
incompressible fluid is the group of all area preserving diffeomorphisms of the fluid container, 
and that solutions of the Euler equations are geodesies on this group with respect to a certain 
kinetic energy metric, characterized by the inner-product J (u- v) dx for two divergence- 
free vector fields u and v. The system (Q) also has this geometric property, but now the 
metric is instead characterized by / (u ■ v + 2ct 2 Def(n) • Def(v)) dx, where Def(u) is the rate 
of deformation tensor (Vm+ (Vm) t )/2.§ Equations @) thus preserve the Hamiltonian 
structure of the Euler equations. In particular, vorticity remains pointwise conserved by 
the smooth Lagrangian flow r]f so that q(r]f (x)) = qo{x), the vorticity momenta / q p dx are 
conserved, and so the Kelvin circulation theorem remains intact! as well. 

Since, in each of the three different scenarios — averaged Euler, vortex-blob method, and 
inviscid second-grade fluid — the essential modification of the original equations is a change 
of the advective nonlinearity of fluid dynamics, we will now consider forced-dissipative sim- 
ulations of the system (|J) to demonstrate the effect of such an inviscid modification. If we 
use a vorticity-stream function formulation, the evolution of both the Euler and averaged 
Euler sytems can be represented by 



where u = Aip, J is the Jacobian, F the forcing, and D the dissipation. (The Euler system 
corresponds to a = 0.) Our numerical scheme consists of a fully dealiased pseudospectral 
spatial discretization and a (nominally) fifth-order, adaptive timestep, embedded Runge- 
Kutta Cash-Karp temporal discretization of (|5]) (see Nadigai for details). With such a 
scheme, among the infinity of inviscid (F — D — 0) conserved quantities for (||), the only 
two conservation properties that survive are those for the kinetic energy E^i and enstrophy 
Z H 2 given respectively by 




(5) 
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(6) 



In the forced-dissipative runs to be considered, the forcing F is achieved by keeping the 
amplitudes of modes with wavenumbers in the small wavenumber band 10 < k < 10.001 
constant in time. The dissipation, D, is a combination of a fourth order hyperviscous 
operator and a large-scale friction term: D = 5ip — (— z/A) 4 u, as has been used in numerous 
previous studies of two-dimensional turbulence. The form and value of the forcing and 
dissipation are held exactly the same for all the runs to be presented, irrespective of the 
resolution and the value of a. 

On the one hand, it could be argued that since the energy and enstrophy that are 
conserved (in an unforced-inviscid setting) are E H i and Z H 2 respectively, it is their dynamics 
which is of primary importance. On the other, it could be argued that in the context of (^|), 
the interest in small scales is only in so much as it affects the larger scales and to that extent 
a has no primary significance and that it is really the large scale components of energy and 
enstrophy in (|6|) that are of primary interest. While both these points of view are reasonable, 
in this short article, we proceed with the latter and concern ourselves with the dynamics of 
the usual kinetic energy and usual enstrophy as given by 



Fig. [I] shows the evolution of the kinetic energy E with time for four different values of k a . 
For these computations, 512 physical grid points were used in each direction, resulting in, 
after accounting for dealiasing, a maximum, circularly-symmetric wavenumber, k max , of 170. 
The four runs correspond to k Q of oo (dissipative Euler) and 42, 21, and 14. This figure 
shows that for identical forcing and dissipation, the tendency with increasing a (equivalently 
decreasing k a ) is to achieve an overall balance which makes the flow less viscous. 

While the kinetic energy of the runs with different a shows a definite trend (increasing 
with increasing a), such is not the case with the enstrophy shown for the same four cases 
in the inset of Fig. |l|. Here interestingly, all the runs with nonzero a seem to display 
approximately the same level of enstrophy which is lower than for k a = 0. This indicates: 

• That the small-scale behavior is quite different when a = and when a is nonzero (as 
noted in previous studies^!'!), but that this difference is not sensitively dependent on 
the value of a for the interesting range of values of a. 

• That the more significant change with a is the behavior of the large scales. 

Therefore, to further examine the nature of this (reduced-viscous) behavior of the large 
scales, we examine the energy-wavenumber spectra in Fig. Here, the average of the 
one dimensional energy spectrum E(k) between times 5 and 20 is plotted against the scalar 
wavenumber k. Figure ^| shows that the reduced- viscous behavior for increasing a is achieved 
by systematically increasing the energy in modes larger in scale than the forcing scale and 
decreasing the energy in modes smaller in scale compared to the forcing scale. 

The larger energy content in the larger scales implies an enhancement of the inverse 
cascade of energy of two-dimensional turbulence by the nonlinear-dispersive modification of 




(7) 
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the advective nonlinearity when a > in (|5|). So also, the decreased energy content in the 
smaller scales is attributable to the same nonlinear-dispersive modification. In the following, 
we give a simple dimensional argument to explain the observed behavior. For this, consider 
the governing equations in the form (||). In close analogy with the classical picture for the 
inertial ranges of two-dimensional dissipative Euler equations,! a Kolmogorov-like cascade 
picture for (|j) shows that the inertial range consists of two subranges, the enstrophy cascade 
subrange where there is a down-scale cascade of the enstrophy defined in (|B|), and the 
energy cascade subrange where there is an up-scale cascade of the E H \ energy defined in 
([]). E H i and Z H 2 are the relevant energy and enstrophies since these are the ones which 
are conserved in an inviscid and unforced case. 

If we assume that the wavenumber k a only appears in the Helmholtz operator, as it does 
in the governing equations, then 

• In the enstrophy cascade subrange, 

E(k) ~ (3 a H2 k\ (8) 



where (3 H 2 is the rate of dissipation of Z H 2 enstrophy, and a and b are exponents to 
be determined by dimensional analysis. If L and T are characteristic length and time 
scales in the enstrophy cascade subrange, (0) implies 



L 3 T -2 = T -3a ft + a 2 L ~2y a £ 



from which a = 2/3. However, even in the enstrophy cascade subrange, the value of 
b depends on the the ratio a/L. For a <C L, of course, b = —3, and the classical! 
E(k) ~ k~~ 3 is recovered; when a ^> L, E(k) ~ k s". Finally, when a is comparable 

17 

to L, it is easy to see that E{k) decays faster than for Euler, but slower than k~~ (as 
may be seen in Fig. 0). 

In the energy cascade subrange, 

E(k)~e a Hl k\ (9) 



where e#i is the rate of dissipation of E H i energy, and a and b are exponents to be 
determined by dimensional analysis. If L and T are characteristic length and time 
scales, now, in the energy cascade subrange, (pi) implies 



L 3 T -2 = T -3a f x + a 2 L -2\ j2a-b ^ 



from which a = 2/3. Again, even in the energy cascade subrange, the value of b 
depends on the the ratio a/L. For a <C L, of course, b = — |, and the classical! 
E{k) ~ k~s is recovered; when q>L, E{k) ~ k~ 3 . When a is comparable to L, it is 
easy to see that the inverse cascade of energy is enhanced, that is E(k) increases with 
decreasing k faster than k~s (Euler) but slower than k~ 3 . 
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While the effects of the enhancement of the inverse cascade of energy is clear in Fig. |2|, 
we defer the verification of the asymptotic values of the exponent b to later studies when 
we can afford much larger simulations with a good dynamic range in each of the inertial 
subranges. 

The steeper fall-off of the energy spectrum with k in the enstrophy cascade range of 
wavenumbers when a > 0, compared to Euler may, at first, suggest that a coarser resolution 
may be sufficient to resolve the flow when a > (for the same forcing and dissipation). 
However, this is not the should be clear from Fig. |3|. In this figure, the spectra of 

the cases previously discussed is replotted together with the corresponding spectra when the 
resolution is reduced by a quarter (k max = 128) and a half (k max = 85). (The spectra for the 
different values of a are offset to improve clarity.) The degree of non-resolution of the flow 
due to the reduced resolution is indicated by the deviation of that spectrum from that for 
the fully resolved case. With a 25% reduction in resolution, the flows are almost resolved 
for all values of a, while with a 50% reduction, the flows are not fully resolved anymore. 
Importantly, the degree of non-resolution is independent of a to the lowest order. 

Besides their use in describing mean motion, the averaged Euler equations have arisen 
independently in at least two other contexts — second grade polymeric fluids and vortex blob 
methods. In this note, we make two observations that are likely to be of fundamental im- 
portance in understanding the relevance of these models in describing more realistic flows: 
While it has been previously noted that with these equations, nonlinear interactions at scales 
small compared to a are suppressed, we have shown here that the modification of the non- 
linear advection term in these equations also leads to an enhancement of the inverse cascade 
of energy in two dimentions — a characteristic feature of two-dimensional turbulence. This 
in turn implies (1) an overall reduced-viscous behavior and (2) a significant modification of 
the dynamics of scales larger than a, both reminiscent of the phenomenon of drag reduction 
in a turbulent flow when a dilute polymer is added (e.g., see [12] and references therein). 
Furthermore, we point out that the limiting of the energy spectrum at small scales due to 
a does not, in itself, allow the flow to be resolved on a coarser grid. 

The latter notwithstanding, we remark that the averaged Euler equations are useful in 
better understanding the limit of inviscid fluid flow, since the averaged Euler equations 
with viscosity, unlike the Euler equations, converge regularly to the solutions of the inviscid 
system.0 That is, for an arbitrary but fixed time interval, we can choose a small enough so 
that the solution of the averaged Euler equations are uniformly within any a priori chosen 
error of the Euler equations^ and then consider the zero viscosity limit of the viscous, 
averaged Euler equations. 

Finally, we note that as a theoretical model of fluid turbulence, the averaged Euler 
equations in two dimensions possess other remarkable features. For example, for any initial 
condition and fixed time interval, one may choose the number of modes k max large enough 
so as to be arbitrarily close to the exact solution of the averaged Euler equations without the 
addition of viscosity. For such large k max , and in simulations of unforced decaying turbulence, 
the averaged Euler equations exhibit a fundametal feature of two-dimensional turbulence: 
a sharp decrease in enstrophy Z during the first few large eddy turnover times. This is 
extremely interesting, because, while it is necessary to add viscosity to the Euler equations 
to obtain similar behavior, the averaged Euler equations can reproduce this behavior while 
exactly conserving an energy. We shall report further on such inviscid simulations in future 
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FIG. 1. The evolution of kinetic energy E with time for k a = oo (solid line), k a = 42 (dotted 
line), k a = 21 (dashed line), and k a = 14 (dot-dashed line). An increase in a, for identical forcing 
and dissipation, results in an overall reduced viscous behavior. In the inset is shown the evolution 
of enstrophy Z for the same time interval and for the same four values of a and with the same 
line types as for kinetic energy. While there is a significant difference between zero and non-zero 
a cases, the dependence on the actual value of a itself is rather weak. 
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FIG. 2. Stationary wavenumber-energy spectra (log-log scale) for the forced-dissipative simu- 
lations of the averaged Euler equations with zero and nonzero a. k a = oo (solid line), k a = 42 
(dotted line), k a = 21 (dashed line), and k a = 14 (dot-dashed line). The insert shows the same plot 
with a linear-linear scale for the first ten wavenumbers. The enhanced inverse cascade of energy 
and the suppressed energy level at smaller scales with increasing a is evident. 



10 




1 10 100 

k 



FIG. 3. The spectra for the four cases in Fig. Q are replotted along with the spectra for the 
same four cases with resolution reduced by 25% and 50% in each direction. The sets of spectra for 
each q are offset by a decade each to improve clarity. The degree of non-resolution of the flow with 
the reduced resolution is indicated by the difference between that spectrum and the spectrum for 
the fully resolved case. With a 25% reduction in resolution, the flows are almost resolved, while 
with a 50% reduction, the flows are not fully resolved anymore. The degree of non-resolution is 
independent of a to the lowest order. 
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